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ABSTRACT 

We investigate the role of the bulk Compton-scattering process in the internal 
shock model. We numerically show that the radiation field created by internal 
shocks does not affect the efficiency of the energy conversion in the shocks even 
if fireball outfiows have highly nonuniform velocity distributions. However, the 
bulk Compton-scattering produces an additional high energy component espe- 
cially when variability in the outfiows arises because of a modulation of the mass 
injected into a constant energy fiow. We show that an isotropic energy of lO^'^"^^ 
ergs could be radiated in the range of lOOMeV-GeV through the scattering pro- 
cess, provided that the shell's Lorentz factor varies between 10 and 10^. 

Subject headings: gamma rays: bursts — shock waves — radiation mechanisms: 
non-thermal 



1. Introduction 

Gamma-ray bursts (GRBs) and afterglows are well described by the relativistic fireball 
model, in which an explosive fiow of relativistic matter is released from a compact central 
engine (e.g. a massive stellar collapse or a compact stellar merger). The kinetic energy 
of the fiow is converted into internal energy by shocks, and it is radiated via synchrotron 
process. These shocks can be due to collision with the ambient medium (external shocks) or 
shocks inside the fiow itself due to nonuniform velocity (internal shocks). The highly variable 
temporal structure observed in most of bursts indicates that internal shocks are responsible 
for the production of GRBs. As internal shocks do not dissipate all the kinetic energy, the 
remaining energy is dissipated later by external shocks and produces an afterglow. 

Since the energy radiated in the afterglow phase is comparable or even smaller than 
that in the prompt gamma-ray phase (Kumar & Piran 2000; Freedman & Waxman 2001; 
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Lloyd- Ronning & Zhang 2004), GRB outflows should be highly irregular. If we represent the 
irregular flow by a succession of relativistic shells with random Lorentz factors (Kobayashi, 
Piran & Sari 1997; Daigne & Mochkovitch 1998; Guetta, Spada & Waxman 2001), the 
distribution of the Lorentz factors should have a large dispersion to achieve a high conversion 
efficiency in the internal shocks (Beloborodov 2000; Kobayashi & Sari 2001). 

Recently, Lazzati, Ghiscllini and Cclotti (1999) pointed out that, if the Lorentz factor 
distribution is very broad, the radiation drag process could affect the conversion efficiency 
and spectrum of the internal shock emission (see also Gruzinov & Meszaros 2000; Ghisellini 
2003; Ghisellini et al. 2000). When the Lorentz factors of shells arc distributed in a large 
interval, fast shells would move in the radiation field produced by the previous collisions, 
and thus Compton-scatter these photons (bulk Compton-scattering) . The fast shells could 
be decelerated before colliding into other shells. This is of significant interest, because high 
energy 7-rays could be produced as a consequence of the internal shock model. Observations 
of high energy radiation will provide valuable information about GRB outflows and their 
central engines. In this Letter, we study the bulk Compton process by using a numerical 
random shell model. We discuss the implication to the Gamma-Ray Large Area Space 
Telescope (GLAST). 



2. Bulk Compton-Scattering 

In the internal shock model, central engines are assumed to produce many relativistic 
expanding shells with a random distribution of Lorentz factors. A collision of two shells is 
the elementary process to dissipate the kinetic energy of the shells (e.g. Kobayashi et al. 
1997). A rapid shell catches up with a slower one and the two merge to form a single one. 
Electrons are heated by shocks, and the internal energy is then radiated via synchrotron 
emission. 

Consider such a merged shell with a Lorentz factor 7, emitting an internal energy E. 
When a fast shell with a Lorentz factor F ^ 7 collides into the merged shell, the fast 
shell scatters the photons emitted from the merged shell. In the comoving frame of the 
merged shell, the internal energy E' ~ E/'j is emitted isotropically. Then, the fast shell is 
moving with the relative Lorentz factor F' ~ r/27 into the photon field, scattering a fraction 
~ min(l, r) of these photons to energies higher by a factor of ~ 2F'^ where r is the optical 
depth of the fast shell. The factor of ~ 2 in the boost factor takes into account that the 
distribution of the photons which the fast shell scatters is not isotropic, and that head-on 
scatterings are favored (GhiseUini et al. 2000). Using the mass of the fast shell M, we can 
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estimate the energy extracted from the fast shell as 




dE. 



(1) 



The fast shell should have a smaller Lorentz factor F^ec when it catches up to the merged 
shell. It is given by 



A fraction of the kinetic energy ~ Mc^(r — Tdec) is radiated away through the scattering 
process. 

In the standard fireball model, the optical depth r of a fireball is close to unity r < 1 
around the internal shock radius, because the parameters are chosen so as to make r < 1 
enabling a non-thermal spectrum, and the "preferred" model uses larger baryon loads (e.g. 
Gruzinov & Mcszaros 2000). Therefore, if merged shells radiate a significant fraction of 
their kinetic energies through the internal shock synchrotron process, as it is assumed in 
the standard model, a large velocity contrast in the fiow F » 7 could lead to an efficient 
production of high energy 7-rays. 



We perform Monte Carlo simulations, and discuss how the bulk Compton-scattering 
affects the evolution of multiple shells. We estimate the conversion efficiency from the kinetic 
energy of shells to high energy radiation in the bulk Compton-scattering. 



We consider a flow consisting of shells. Each shell is characterized by four variables: 
a Lorentz factor ji, a mass rrii, a width /j and a position Ri. We assume that the initial 
Lorentz factor of each shell is distributed uniformly in logarithmic space between log 'jrnin 
and log jrnax- The initial masses are assumed to be correlated with the Lorentz factors as 
m oc 7*. For s = — 1 and 0, the shells initially have equal energy and equal mass respectively, 
and for s = 1 the shells initially have equal density under an assumption of the equal shell 
width. We assume a constant value / for the initial widths and the initial separations between 
the shells. Then, the initial positions of the shells are — 2{i — 1)1. 




(2) 



3. 



Numerical Shell Model 



3.1. 



Outline of the Model 



The evolution of the shell system in time is basically the same as in Kobayashi et 
al (1997), but the radiation drag effect is taken into account. Because of the relativistic 



4 



beaming effect, tfie radiation from a merged shell with a Lorentz factor 7 is concentrated 
in the forward direction. A large fraction of photons are emitted around 6 ^ l/'j, and they 
propagate to the radial direction with a velocity of ~ c - cos(l/7) or equivalently an effective 
Lorentz factor of ~ 7. We regard the radiation from each collision as a photon packet with 
the Lorentz factor of the merged shell 7. Since merged shells will further collide with other 
shells, the Lorentz factors of the mergers would significantly change during the evolution. 
The Lorentz factors of photon packets are assumed to be constant after the productions. 

Wc turn now to describe the evolution of the multiple shell system. There will be 
numerous collisions between the shells. We denote by the index j the j-th collision which 
takes place at a time tj. Given the velocities of the shells (3i and their separations at we 
calculate for all pairs {i, i + 1) the collision times. We then find the minimal collision time 
5tj. After rearranging shells from earlier positions Ri{tj-i) to Ri{tj) = Ri{tj-i) +c(3i5tj, we 
consider the radiation drag effect. Each photon packet is characterized by a position r„, an 
energy En and a Lorentz factor rjn {n = 1, .., j — 1). To estimate the radiation drag effect 
on shell i, we find all the packets {k} which the i shell overtakes between tj-i and tj. The 
Lorentz factor of the i shell at tj is evaluated by 



where Tj is the optical depth of the i shell at tj. and wc take a sum over the subset {k}. The 
difference of the kinetic energies rriiC^ [7j(tj_i) — 7i(tj)] is emitted from the i shell as high 
energy photons. Since a fraction of the photons in the relevant packets {k} are scattered by 
the i shell, we reduce the energy as E^itj) = [1 — min(l,ri)] E]^{tj_i). 

We estimate the internal energy produced by the j'-th collision by using the energy- 
momentum conservation (eq. (2) in Kobayashi et al. 1997). We assign 1/3 of the internal 
energy to the j-th photon packet produced in this collision if the optical depth of the merged 
shell is less than unity. Wc assume that the remaining internal energy is transformed to 
the kinetic energy of the merged shell. If the optical depth is larger than unity, the whole 
internal energy is converted to the kinetic energy of the merger. We find the next collision 
and proceed until there are no more collisions, i.e. until all the shells have merged to form 
a single shell or until the shells are ordered with increasing value of the Lorentz factors (the 
outermost and innermost are fastest and slowest, respectively). 

Photon packets can boost shells in principle when they catch up with slow shells. How- 
ever, the boost effect is of the order of rdE instead of {V/^YrdE which is the order of the 
radiation drag effect (see eq. (1)). We find that the boost effect is negligible in the evolution 
of the multiple shell system. In this letter, we discuss only the radiation drag effect. 




(3) 
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3.2. Efficiency 

For N = 300 shells, jmin = 10, l/c = 1 sec and the initial kinetic energy Ekin,iso = 
'^miC^'-fi = 10^^ erg, we have evaluated the conversion efficiencies of internal shocks and 
the bulk Compton-scattering. The mean efficiencies and their standard deviations for 100 
reahzations are hsted in Table 1. In the third column, "No Scattering" indicates the internal 
shock efficiencies obtained without the scattering process. We can see that the efficiencies 
of the internal shocks are not influenced by the radiation drag effect significantly even when 
the velocity dispersion is large ^max/lmin ^ 1- 

While the internal shock process is most efficient in the equal mass cases, the bulk 
Compton-scattering is most efficient in the equal energy cases. Recently Nakar and Piran 
(2002) studied time scales in GRB temporal profiles. Their results imply that the variability 
in wind's Lorentz factors arises because of a modulation of the mass injected into a constant 
energy flow (equal energy shells). If this is the case, a significant amount of energy, corre- 
sponding to 1 — 10% of the prompt gamma-ray energy (GRB energy), might be radiated as 
high energy 7-rays through the bulk Compton-scattering. 

A large number of collisions ~ happen during the evolution of the multiple shells. 
The ffist collisions occur between the slowest shells ~ 7mm and fast shells ^ 7min around 
R ~ ImiJ- equal energy cases, the Lorentz factors of the merged shells produced by the 
ffist collisions are comparable to that of the slowest shells ~ 7r„m (eq. (1) in Kobayashi et al. 
1997). This property allows other fast shells to catch up with the merged shells at small radii 
~ 27^j^/. Since the fast shells still have larger optical depth oc R"'^ at the small collision 
radii, and the relative Lorentz factors are larger, the photons emitted from the first mergers 
are efficiently scattered into high energies. In equal density cases, the Lorentz factors of the 
first mergers are close to ~ ^max- It is difficult for other shells to catch up to the mergers. 
Since such collisions happen at large radii and the relative velocities are small, the bulk 
Compton-scattering is expected to be inefficient. In equal mass cases, the Lorentz factors of 
the first mergers take intermediate values ~ {'^minlmaxY^'^ ■ The efficiency of the scattering 
process is also expected to be intermediate. 

We classify shell collisions occurring in the optically thin region into four cases: (I) 
Collisions happen between two shells, both of which have never experienced collisions. Since 
we are interested in the radiation drag effect, we count only collisions in the optically thin 
region. Even if a shell has collided with another shell in the optically thick region, we regard 
the merger as an unexperienced shell. (II) While the outer slow shell went through collisions 
at least once, this is the first coUision for the inner rapid shell. (Ill) The outer slow shell has 
never interact with other shells, and the inner rapid shell has experienced collisions at least 
once. (IV) The other cases. Assuming the same parameters as in Table 1 and ^max — 10^, we 
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have performed numerical simulations for the three mass distributions, and we have analyzed 
all shell collisions in the optically thin region. We summarize the results for 100 realizations 
in Table 2. The radiation drag effect is expected to be most significant in type II collisions. 
Actually, in the equal energy case for which the bulk Compton-scattering is most efficient, 
about 50% of collisions are classified into the type II. On the other hand, the fraction of type 
III collisions is larger than that of the type II in the equal density case. These results are 
consistent with the above discussion. 

Efficient bulk Compton-scattering requires the optical depth of shells to be close to unity 
r < 1 at internal shock radii. In fig 1, we plot efficiencies as functions of the variability time 
scale //c, and we discuss how the optical depth is sensitive to the choice of the variabiUty 
time scale. If it is very short Z/c ~ 1ms, most collisions happen at small radii. Since the 
shells still have large optical depth, the internal shock emission is highly suppressed, and 
consequently only a small amount of energy is up-scatted to high energies. As a larger value 
of //c is assumed, a larger fraction of collisions happen in the optically thin region. AroTind 
i/c ~ a few sec, all the collisions happen in the optically thin region, and the internal shock 
efficiency (thin solid line) is saturated. The efficiency of the bulk Compton-scattering (thick 
solid line) also increases as a larger //c is assumed, because the internal shocks create thicker 
photon field. However, if the variability time is very long > 10s, the optical depth of the 
shells at collision radii is too small to up-scatter the photon field effectively. For our typical 
parameters, the efficiency of the bulk Compton-scattering peaks around Z/c ~ 0.1 — 1 sec, 
which are comparable to the typical variability time scales observed in GRB profiles. 

All the numerical results obtained so far are for Ekin,iso = 10^"^ ergs, which is the average 
value of what we derive from observations. The distribution of isotropic GRB energies spans 
at least two orders of magnitude. Evaluating the efficiencies for 10^^ ergs (dashed fines in fig 
1) and for 10^^ ergs (dashed-dotted lines), we found that the asymptotic values of internal 
shock efficiencies ~ 11% and the peak values of scattering efficiencies ~ 3% are the same, and 
that the efficiency curves just shift to left or right. As we can expect from the optical depth 
of shells r oc Ekin,iso/i'^, the saturation time scale for lO^^ergs, at which the internal shock 
efficiency is saturated and the scattering efficiency peaks, is one order larger than that for 
lO^^ergs. For the same distribution of Lorentz factors, the efficiency of internal shocks does 
not depend on Ekin,iso ff collisions happen in the optical thin region. The deceleration 
effect dT by radiation field is larger for a smaller baryon loading (a smaller Ekin^iso)- However, 
the energy extracted from shells depends on the mass of shells (baryon loading) also, and it 
is determined by the photon field energy as Mc^dT oc dE. The peak values of the scattering 
efficiencies should have similar values, because the internal shock efficiencies are the same at 
the saturation time scales. 
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We have assumed that 1/3 of the internal energy produced by each coUision goes into 
radiation. Assuming a smaller value of 0.1 or 0.01, we re-evaluated the efficiencies in the 
equal energy case of ^{max = 10^ in table 1. The internal shock efficiencies arc 3.6 ± 0.2% 
and 0.37±0.02% for the conversion factors of 0.1 and 0.01, respectively, while the scattering 
efficiencies are respectively 1.1 ± 0.3% and 0.14 ± 0.04%. If less photons are produced, less 
photons would be scattered. The ratio of internal shock and scattering efficiencies is almost 
constant. 

We found that the internal shock process is more efficient than the bulk Compton- 
scattering. The main reason is that internal shocks occur at a rather wide range of radii. 
Since the optical depth of shells rapidly decreases as r oc R^^, the Compton-scattering pro- 
cess is effective only at small radii at which only a small fraction of the kinetic energy of 
the flow is converted to photon fleld. We have analyzed collisions occurring in an equal 
energy case with the same parameters as in fig 1 and l/c ~ 1.5s (the scattering process is 
most efficient for this /) , and we obtained the following results. The internal shock emission 
starts around i? ~ 4 x lO^^cm, and the shocks at i? < a few lO^^cm contribute effectively 
to the total internal shock radiation of ~ lO^^ergs. 10%, 50% and 80% of the radiation 
are emitted within R ~ 10^^ cm , 5 x lO^'^cm and 2 x lO^^cm, respectively. The mini- 
mum internal shock radius ~ 4 x 10^'^cm is consistent with the photosphere radius Rph 
at which mergers becomes optically thin. It is given by Rph ~ {o-t Ekin,iso/ ^^n N'jrnin'nT'pC^y^'^ 
~ 3.4x 10^^(£;jfci„,iso/10^^ergs)^/2(^iV/300)-^/2(^7^in/10)-^/2 cm. The bulk Compton-scattering 
process starts around ~ 2i?p^, at which the optical depth of scatter shells with F » 7^j„ is 
already much smaller than unity, because their masses are smaller than the masses of the 
mergers. Consequently, the Comptonization parameter y = rx (relative Lorentz factor)^ is 
on the order of unity at i? < lO^^cm and it decreases rapidly at larger radii. 10%, 50% and 
80% of the bulk Compton emission are produced within i? ~ 9 x lO^^cm, 2 x 10^"^ cm and 
4 X lO-'^^cm, respectively. Although the scattering process is very effective y ~ 1 at small 
radii R < 10^"^cm, the energy of the photon field created by the internal shocks at the small 
radii is ~ 10% of the total internal shock radiation, and then the efficiency from the kinetic 
energy to high energy radiation in the bulk Compton-scattering is of the order of 1%. If we 
assume even larger contrasts of Lorentz factors (7^02: = 10^ or larger), scattering shells have 
larger Comptonization parameters ?/ ^ 1 at the first scattering radius ~ '2Rph- However, 
the scattering emission with large relative Lorentz factors is suppressed by the Klein- Nishina 
effect. The resulting efficiency of the scattering process is about 4-5 % for ^max — 10^~^. 

The spectra of the high energy radiation produced by the bulk Compton-scattering are 
shown in fig 2 for two example cases. Assuming that all the primary photons produced by 
shell collisions have the typical GRB frequency of 100 keV, we have integrated up-scattered 
photons to construct the spectra. Bursts of high energy emission are expected to occur at the 
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same time with the prompt 7-ray radiation (GRBs). The durations should be comparable 
to these of the associate GRBs. If the Lorentz factors of shells are distributed between 10 
and 10^, the typical frequency of the bulk Compton emission is about lOOMeV-GeV. Since 
the primary photons are scattered to frequencies higher by a boost factor oc (r/7)^ where 
r and 7 are respectively the Lorentz factors of the scatter and emitter shells, the typical 
frequency of the bulk Compton emission for jmax = 10^ is lower, and the emission peaks 
around 10-100 MeV. 



4. Discussion 

We have studied the bulk Compton-scattering in the internal shock model by using 
Monte Carlo simulations. We have shown that the radiation drag does not affect the conver- 
sion efficiency in the internal shocks even if fireballs have a very large velocity irregularity. 
The bulk Compton-scattering is most efficient when an inner engine produces shells with 
comparable energy but very different Lorentz factors. Such equal energy cases are favored in 
a recent study on GRB temporal profiles (Nakar & Piran 2002). We have obtained about a 
few % conversion efficiency of the bulk Compton-scattering, provided that the shell's Lorentz 
factor varies between 10 and 10^. This implies that Eiso ~ lo^'^"^^ ergs might be radiated as 
lOOMeV-GeV 7-rays. 

The next generation Gamma- Ray Large Area Space Telescope (GLAST), scheduled for 
launch in 2007, will provide unprecedented sensitivity in the energy range of about 20 MeV 
to about 300 GeV. GLAST should be able to detect the high energy emission produced by 
the bulk Compton-scattering if it is associated with GRBs. GLAST observations will provide 
a variable clue to understand the nature of GRB outflows and the central engines. 

Although the synchrotron self-inverse Compton-scattering also can produce high energy 
7-rays, the typical frequencies (> GeV) could be much higher than that of the bulk Compton 
emission especially when fireballs have large velocity fluctuations as we have assumed in this 
Letter. Recent studies suggest that flreballs could be highly magnetized (Zhang et al. 2003; 
Kumar & Panaitescu 2003; Fan et al. 2002). If this is the case, the self- inverse Compton 
emission could be suppressed. However, the bulk Compton-scattering could still produce 
high energy 7-rays as long as fireballs are not dominated by Poynting-flux. 

This work is supported by the Eberly Research Funds of Penn State, by the Center for 
Gravitational Wave Physics, funded under cooperative agreement by NSF PHY 01-14375, 
by a NASA Swift Cycle 1 GI program. 
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Table 1. Efficiency 


Mass Distribution 




Internal Shocks (No Scattering) 


Bulk Compton 


Equal Energy 


10^ 


8.1±0.4 % (8.2 ± 0.4 %) 


0.7 ± 0.2 % 


Equal Energy 


10^ 


10.5±0.4 % (10.9 ± 0.4 %) 


2.7 ±0.6 % 


Equal Mass 


10^ 


24.0±1.4 % (24.0 ± 1.4 %) 


(5.8 ±3.7) X 10-2 % 


Equal Density 


10^ 


14.5±1.4 % (14.5 ± 1.4 %) 


(1.0 ±0.9) X 10"^ % 



Table 2. The basic types of two shell coUisions 



Mass Distril")ution 


Ty])e I 


Type II 


Type III 


Type IV 


Equal Energy 
Equal Mass 
Equal Density 


23.7±1.4 % 
32.1±1.3 % 
33.8±1.4 % 


50.4±2.8 % 
24.4±2.3 % 
16.5±2.2 % 


4.0±1.2 % 
13.3±2.1 % 
17.9±2.5 % 


21.9±1.5 % 
30.2±1.4 % 
31.9±1.5 % 
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Fig. 1. — Efficiency vs. variability time scale. Bulk Compton efficiency (thick) and internal 
shock efficiency (thin). Ekin,iso = 10^^ ergs (dashed), 10^^ ergs (solid) and 10^'' ergs (dashed- 
dotted). 300 equal energy shells, ■jmin — 10 and jmax ~ 10^. 
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Fig. 2. — Spectra of bulk Compton emission, ^max — 10^ (dashed line) or ^max — 10^ (solid 
line). The other parameters are Ekin,iso = 10^^ ergs, 300 equal energy shells, 7min = 10 and 
l/c — 1 sec. The cutoff of the solid line around a few GeV is due to the Klein- Nishina effect. 



